Cluster virial expansion for nuclear matter within a quasiparticle statistical approach 
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Correlations in interacting many-particle systems can lead to the formation of clusters, in partic- 
ular bound states and resonances. Systematic quantum statistical approaches allow to combine the 
nuclear statistical equilibrium description (law of mass action) with mean-field concepts. A chemical 
picture, which treats the clusters as distinct entities, serves as an intuitive concept to treat the low- 
density limit. Within a generalized Beth-Uhlenbeck approach, the quasiparticle virial expansion is 
extended to include arbitrary clusters, where special attention must be paid to avoid inconsistencies 
such as double counting. Correlations are suppressed with increasing density due to Pauli blocking. 
The contribution of the continuum to the virial coefficients can be reduced by considering clusters 
explicitly and introducing quasiparticle energies. The cluster-virial expansion for nuclear matter 
joins known benchmarks at low densities with those near saturation density. 
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I. INTRODUCTION 

Recently considerable efforts have been made to elaborate the nuclear matter equation of state (EoS) in a wide 
range of baryon density n, temperature T, and proton fraction Y p (or neutron-proton asymmetry 6=1 — 2Y p ) [TH2]. 
We consider in this work warm dilute matter (T < 20 MeV, n < uq with no ~ 0.16 fm~ 3 being the nuclear saturation 
density) . In this region the simple model of an ideal Fermi gas of particles has to be improved by including correlations 
between the particles. It has to be emphasized that the notion of correlations depends on the reference state. For 
an ideal Fermi gas, the many-body state is given by a Slater determinant of single-particle plane waves without any 
correlations beyond those originating from quantum statistics. Explicit correlations become less significant near the 
saturation density because correlations are increasingly blocked with increasing density and Fermi energy and the 
many-body state can be considered as a system of uncorrelated quasiparticles. At low densities an explicit treatment 
of correlations is essential. In particular, clusters can dominate the composition of matter at finite density at low 
temperatures. Different approximations have been considered to treat the formation of clusters in the low-density 
limit of the nucleon gas, such as the nuclear statistical equilibrium (NSE) approach 6, 8J or virial expansions [9l412J . 
Here we want to show that a consistent treatment of correlations (bound as well as scattering states) can be given 
within a quantum statistical approach. 

The main question is to match quasiparticle concepts with cluster concepts, as e.g. embodied in the law of mass 
action at low densities. This has been discussed for the two-nucleon problem within a generalized Beth-Uhlenbeck 
approach [T3] which was developed originally for electron-hole plasmas [T3]. A general expression for the second 
virial coefficient was obtained that contains in addition to the contribution of bound states also those of scattering 
states. As an important result, it has been demonstrated that the contribution of the continuum of scattering states 
is modified when quasiparticles are introduced. Larger clusters consisting of A nucleons were introduced via a cluster 
decomposition of the nucleon self energy [TBI . A quasiparticle treatment of nuclei was considered via an in-medium 
Schrodinger equation, where, in addition to the single-nucleon self-energy shift, also the Pauli blocking was considered. 
In-medium quasiparticle energies for the ground states of light nuclei A < 4 have been given recently |16j as function 
of temperature T, baryon density n, proton fraction Y p , and the center-of-mass momentum P. The solution of the 
in-medium few-body problem would also give excited states as well as the scattering states. 

Quantum statistical approaches are based on perturbation expansions. Within the Green function method, Feynman 
diagrams are introduced, and partial summations are performed. We focus on a special prescription to select relevant 
contributions of the perturbation expansion that describe correctly the formation of clusters in the low-density limit. 
The concept to include bound states on the same footing as new particle species that can react is called the chemical 
picture. On the other hand, a fundamental quantum statistical approach of interacting elementary particles is the 
physical picture. From this the chemical picture is obtained if in addition to diagrams with single-particle (quasi- 
particle) propagators also diagrams are considered where the single-particle propagators are replaced by ladder-type 
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diagrams that describe the propagation of the A particle cluster. If only the bound state contributions are included, 
one obtains the NSE in the low-density limit. 

If scattering contributions are included in the two-body channel, one obtains the virial expansion. A generalization 
of the virial expansion is also possible within the chemical picture if in any channel of the cluster-cluster interaction not 
only the formation of bound states, but also the contribution of scattering states to the thermodynamical properties of 
the system is considered. Indeed, there are regions in the density-temperature plane where clusters are dominant and 
should be considered as new constituent particles within the chemical picture. In this case, a cluster-virial expansion 
should be performed to include the effects of the continuum. Such an approach is evident from the empirical point 
of view, but not easily derived from first principles. The use of a cluster mean-field (CMF), as well as nearly bound 
states from the continuum correlations can lead to double counting, and these contributions have to be extracted from 
the continuum contribution. 

Recently, such a cluster-virial expansion has been applied to the proton-neutron-a system |12j and was extended 
to include additional light clusters [T7]- Bound states, such as the deuteron or resonances like 8 Be, were treated 
differently, sometimes as part of the generalized second virial coefficient, sometimes as new constituent particles. A 
systematic derivation of a cluster virial expansion that includes as a limit also the NSE has to be performed on the 
basis of a fundamental quantum statistical approach. Contributions occurring in higher-order virial coefficients that 
correspond to bound states or resonances have to be separated since they are explicitly accounted for in the NSE 
including all nuclei in ground and excited states. 

The chemical picture can also be used to derive a cluster-mean field approximation that treats the mean-field effects 
of a correlated medium [TS]. It is expected that the incorporation of mean-field effect into the cluster quasiparticle 
states will also change the contribution of the scattering states in the cluster-virial expansion. This can be understood 
from a fundamental point of view by considering the spectral function in the respective A-nucleon channel. Cluster- 
quasiparticles should be introduced representing the peaks of the A-nucleon spectral function. Then the explicit 
treatment of the contribution of these peaks as quasiparticles will account for a significant part of the total contribution. 
This is clearly seen in the low-density region where a law of mass action can be introduced. The occurrence of bound 
states in the EoS according to the law of mass action is a signature that significant contributions have to be extracted 
from the spectral function of the elementary nucleons. The ordinary quasiparticle picture where the nucleonic spectral 
function is assumed to be sharply (5-like) peaked is no longer justified. 

The inclusion of scattering states is of importance for the equation of state as well as for further properties of nuclear 
systems. Recently, the composition of low-density nuclear matter was investigated [J and it was seen in the quantum 
statistical approach that the contribution of scattering states led to a reduction of the deuteron mass fraction. This 
has to be taken into account in particular at high temperatures. We will give some details here. The main result is 
that the contribution of the continuum can be reduced if mean-field and cluster contributions are already extracted. 

Some results of the chemical picture are given in Sec. [ITJ however a comprehensive treatment cannot be given here. 
We discuss the different effects and propose some approximations. The main aim of this paper is to obtain a more 
systematic treatment in particular of the cluster-virial expansion and to avoid inconsistencies such as double counting. 
As example, the second virial coefficient is considered in Sec. |III| that gives the leading contribution in the low density 
limit. In particular, we discuss the ambiguities connected with the introduction of the bound state contribution of the 
virial coefficient and the consistent introduction of the quasiparticle picture. The relation to the generating functional 



approach that allows for a systematic quantum statistical approach to thermodynamic properties is outlined in Sec. IV 
Conclusions are drawn in Sec. |Vj and details are given in the Appendices. 



II. THE CHEMICAL PICTURE AND CORRELATIONS 



We explain first some general results for many-particle systems. In the chemical picture, we start with a mixture 
of different constituents that can be the elementary particles (atoms, here nucleons) as well as the bound states 
(molecules, here nuclei). The main issue is to design approximations where "elementary" particles and "composite" 
particles are treated on the same footing. This concept has to be incorporated into a quantum statistical approach 
using Green functions techniques. The cluster Green function approach is briefly summarized in App. [XJ 



A. Nuclear statistical equilibrium 

At low densities we can neglect the interaction with exception of collisions where the constituents come close together 
resulting in reactions that establish the chemical equilibrium. These reactions include excitations and ionization, e.g., 
in a plasma. The same considerations apply also for nuclear systems where various nuclei occur that can react. Thus, 
we have as approximation an ideal mixture of different components. 



3 



As a result, the law of mass action is found, with the total particle (baryon) number density (note that astrophysical 
/^-equilibrium is not considered here) 

n(T,fi p ,fi n ) = J2 ^E/a^z^Wz) (1) 

A,Z,v p 

where A, Z denote mass and charge number of a nucleus, respectively, and v indicates the internal quantum state of 
the nucleus. P is the center of mass momentum, Q is the volume of the system, and 

U{E > M) = cxp[(£- M )/T]-(-l)A ( 2 ) 

is the Bose or Fermi distribution function for even or odd A, respectively. Note that at low temperatures Bose-Einstein 
condensation may occur, which is neglected here. We determine the chemical potentials fiA,z = Zfi p + (A — Z)fi n 
only with respect to the kinetic energy of species {A, Z, is}, the binding energy Ba,z,v is considered explicitly. 
In the low-density limit, the energies 

E a!zA P ) = ~ B A,z,v + P 2 /(2m A ^ u ) (3) 

are given by the binding energies Ba,z,v of the isolated nuclei in vacuum with masses itia.z,l> — Zm p + (A — Z)m n — 
Ba,z,v, denoted by the index (0). In the NSE the summation over v concerns only bound states, the contribution of 
the continuum is neglected. Neglecting medium corrections, the ordinary law of mass action will increase the mass 
fraction of bound states like d and even more the a particle with respect to the free proton and neutron fractions 
when the total nucleon number density n increases. This behavior contradicts the expectation that at high densities 
a nucleonic quasiparticle picture is appropriate without the explicit occurrence of clusters. 

For a more fundamental approach that is not limited to the low-density limit, one should use the physical picture 
where some of the constituents are elementary, while others are composite particles. The composite particles are 
obtained as bound states of the elementary particles. The same interaction potential, that leads to the formation 
of bound states in the solution of the many-body Schrodinger equation, determines also the interaction between the 
constituents of the system. (Note that it is common to formulate statistical physics with 'elementary' particles that 
on their part are composed of more elementary particles.) The simple law of mass action with clusters where the 
interaction between the components is neglected (with exception of reactive collisions) works well in the low-density 
region and low temperatures, in contrast to a picture of an ideal mixture of the elementary particles such as the ideal 
Fermi gas of protons or neutrons. 



B. Cluster-virial approach 



We discuss three issues that are of relevance to improve the NSE. The first one refers to the summation over excited 
states that, formally, can also be considered as new, independent constituents. Contributions of the continuum are 
neglected within the chemical picture, and only bound clusters are considered as new 'components'. The limits of this 
simple chemical picture are evident for sharp resonances above the continuum threshold. Long living states should 
be included in a more general treatment if the life time (inverse of the width) is sufficiently long. We will consider 
this problem below within the physical picture using a quantum statistical approach. To obtain the complete second 
virial coefficient, however, the contribution of scattering states has to be included according to the Beth-Uhlenbeck 
approach. 

The second issue refers to the treatment of the interaction between the different components. The chemical pic- 
ture takes interactions into account only to establish chemical equilibrium. Effects of scattering correlations on the 
thermodynamical quantities are neglected. To obtain the full second virial coefficient, the contribution of scattering 
states has to be included. Nevertheless, the chemical picture describes well the low-density, low-temperature region. 
It can be improved taking so-called excess terms in the chemical potential and other thermodynamic variables into 
account. As an example, a simple way to include non-ideal effects is the excluded volume concept 6j. More system- 
atic approaches consider scattering phase shifts due to the interaction between the constituents, as considered in the 
Beth-Uhlenbeck formula. Then, the total baryon number density 

n(T, n P , (i n ) = m(T, fip, fin) + n 2 (T,fi p ,fi n ) + n 3 (T, fi p , fi n ) + ... (4) 

contains a contribution n\ from the individual constituents, i.e. nucleons and bound states of nuclei. This term is 
identical to the NSE result Eq.ffil. The contributions n%, n 3 , ... account for the two-body, three-body, . . . correlation 
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effects in the continuum. In the original Beth-Uhlenbeck formulation there are contributions from bound and scattering 
two-body correlations appearing in n.2. However, in order not to count contributions twice as original constituents in 
Tlx and as bound states in rt 2 , the latter term should contain only the scattering part. It can be expressed through 
integrals with scattering phase shifts 8 C in all channels c for the scattering of nuclei {A, Z, v\ and {A', Z', v'} . Thus 
it is given by 

n 2 (T, Mp)Mn )= £ J2 ^EE ^'V" r dEf A+A ,(E^(P) + E^ A , z+ ^, z ,)^ (5) 

\ /..v V. p 

with energies 



E<® (P) = -B A , ZtV - B A , tZ , y + , (6) 

2{m A , ZtV + rriA>,z',v>) 

and degeneracy factors g c . 

The third issue is related to medium effects. The energies E A 0) Zv {P) and E [ c\P) in the NSE and the virial 

expansion, respectively, will be replaced by quasiparticle energies E A>ZtU (P;T, fi p , fi n ) of the clusters leading to the 
generalized cluster Beth-Uhlenbeck approach that is discussed in the next subsection. These quasiparticle energies 
depend on density and temperature in many-particle systems of finite density. Medium effects become operative if 
the density exceeds about 10 -4 fm -3 . In particular, Pauli blocking will dissolve the clusters at densities about 10~ 2 
fm" 3 [16]. 

The quasiparticle approximation has been used in several recent approaches to calculate the properties of nuclear 
matter [HO [7]. In particular, it is very effective near saturation density but fails in the low-density region where 
clusters eventually become of importance. As a first step, we discuss the inclusion of two-particle correlations and 
the generalized Beth-Uhlenbeck formula, i.e. the n,2 term in Eq.Q, in the following Section. It reproduces the 
correct second virial coefficient but gives also the correct high-density limit. This allows to reproduce the NSE in the 
low-density limit but also the transition to the correct behavior in the high-density region. 

It will be shown, that both effects, the consideration of the interaction between clusters, i.e. of cluster mean- 
field effects, as well as the taking into account of excited states and resonances, will reduce the contribution of the 
continuum to the equation of state, so that appropriate approximations can be obtained. 



C. Generalized cluster Beth-Uhlenbeck approach 



Now we discuss the inclusion of medium effects. The effective degrees of freedom are now quasiparticles with 
self-energies that replace the original constituents. The full evaluation of the imaginary part of the single-nucleon 
self-energies Im Si in the so-called T 2 G\ approximation, see Ref. [13], leads to a generalized Beth-Uhlenbeck equation 
and corresponding EoS. We can extend this approach, originally formulated with only nucleons as basic constituents, 
heuristically to consider nucleon and cluster quasiparticles. A systematic treatment based on a generating functional 



approach is sketched in Section IV 



Considering two-body correlations of nucleons and clusters at most, the total baryon number density 

n(T, fi p , fi n ) = n\ u (T, fi p , fi n ) + n^ u (T, fi p , fi n ) (7) 

receives contributions from single quasiparticles (n, p, d, t, 3 He, 4 He, etc.) and correlated two-quasiparticle contin- 
uum states. The one-body term 

"i U ( T >Mp>Mn) = XI ^ X fA(E A , z ,u(P;T,fi p ,fi n ),fi AtZtl/ ) (8) 

A,Z,v 



P>P1 



Mott 



resembles the corresponding contribution ([!]) in the NSE and ni(T, fi p , fi n ) in the conventional cluster-virial approach 
discussed in subsections |II A| and |II B| However, there are distinct differences. The quasiparticle energies 



E AtZiU (P;T,fi p ,fi n )=E^ z (P)+AE A %(P ] T,fi p ,fi n ) (9) 



depend not only on the cm. momentum P of the nucleon or cluster with respect to the medium but also on the 
temperature and chemical potentials (or equivalently densities) through the medium dependent self-energy shift 
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AE A E Z >V (P\ T, (Ap, fi n ). See App. [Bjand Ref. p5J for a derivation of these shifts. In addition, bound states of clusters 
exist only for cm. momenta P that are larger than the (quasiparticle and medium dependent) Mott momentum Pmom- 
This is mainly a consequence of the Pauli principle that suppresses the formation of clusters due the population of low- 
momentum states by nucleons of the medium. The Mott momentum PMott{T, n n , n p ) indicates the critical momentum 
where the bound state merges with the continuum of scattering states. The Mott momentum becomes larger than 
zero if the density exceeds a critical value, i.e. above this density bound states can exist only for P > PMott(T, n n , n p ) 
[15] , Obviously, Pmou — for nucleons. 

The two-quasiparticle scattering contribution 

qu/rp x A + A' y—^ V"^ * + $A,Z,u;A',Z',u' /, n x 

n 2 {T, Mpj Mn) = 2^ 2^ —q— 1^1^ 9c ^ (10) 

A,Z,v A' ,Z' ,v' S c 



X 



f°° / \ dS 

J dE Ja+A' [E c (P] T, (ji p , fi n ) + E, n A , z + Ha>,z>) 2 sin 2 (S c ) 



is also modified in comparison to the virial result ^ since the energy 

P 2 

E C {P; T, ftp, ii n ) = -B AZiV - B A , !Z ,y + ■ + AE* E (P; T, fi p , M „) (11) 

2(mA,z,i/ + mA>,z',u>) 

contains the medium-dependent shift AE^ E (P; T, (i p , // n ) of the continuum edge as determined by the the self-energy 
shifts of the free constituents. Moreover, there is an additional 2 [sin (S c )} 2 factor that reduces the two-body scattering 
contribution because a part of the two-body correlation effect is shifted to the self-energies of the quasiparticles. 



Note that the two-particle contribution (10) describes binary scattering processes A, A'. A possible bound state 



contribution in that channel is excluded as contribution to the two-quasiparticle scattering contribution (10). The 
contribution of a possible A + A' bound state is already taken into account in the one-body term ([8]). Thus, double 
counting is avoided. 



III. COMPARISON OF APPROACHES AT LOW DENSITIES 



In the low-density limit, medium effects can be neglected. Only neutrons, protons and deuterons are the relevant 
constituents. The standard Beth-Uhlenbeck formula for the second virial coefficient and the virial expansion [9l [10] 
are exact results. They are derived by expanding the thermodynamic functions with respect to the fugacities. It is 
instructive to show explicitly the equivalence with the generalized cluster Beth-Uhlenbeck approach that introduces 
the concept of quasiparticles to account for part of the correlations by introducing self-energies. 



A. Cluster virial approach 



Let us first consider the cluster virial method applying a fugacity expansion up to second order in the chemical 
potentials of protons and neutrons. With the continuum approximation (l/fi)X}p "~ * / d 3 P/(2w) 3 the one-body 
contribution assumes the form 



(p) 

ni = n{ 



An) 



1 ' 1 

with the single nucleon contributions (i — p, n, including degeneracy effects) 

'2[ii 



,(i) 



2 

A? 



exp 



m 



exp 



T 



and the single dcutcron contribution 



,('') 



exp 



+ Mn + B2.1 
T 



(12) 



(13) 



(14) 



A, = y/2ir/(m i T) are the thermal wavelengths of the particles i — p,n,d and B 2 ,\ > is the binding energy of 
the deuteron ground state (identical for all three substates with v = J z = —1,0,1). The two-body term n 2 in 
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this approximation is limited to the two-nucleon scattering contributions. Hence, no nucleon-deuteron or deuteron- 
deuteron correlations are included. In the following, we consider only s-wave contributions for simplicity and have 



exp 



2/Xp 
T 



T 



.3/2.3/2 



exp 



T 



(15) 



with the (continuum) virial coefficients (assuming for strong interactions the symmetry b p 



and b r . 



b - 2 3/2 
= 2 1 ' 2 



dE 




' E 




exp | 


, T 


dE 




' E 


7T 


exp | 





E\ d8\ 



(nn) 



So 



dE 



(pn) 



So 



dS 



(pn) 
3 S t 



dE 



dE 



(16) 
(17) 



with phase shifts in the isospin triplet (T — 1, 1 Sq) and singlet (T = 0, 3 Si) channels. The phase shifts can be taken 
from experiments such that no further model parameters are needed |12) . Instead one can also use for their calculation 
a nucleon-nucleon interaction potential which is adjusted to describe the empirical scattering phase shifts. 



B. Ambiguity of bound state contributions and physical picture 



A partial integration gives the alternative expressions 



bn,?} 



2 3/2 


r°° dE 




' E 


~Y~ 


JO 7T 


exp | 




2 l/2 


r°° dE 




' E 


~~T~ 


Jo ^ 


exp | 





6^(E) + 3S^(E) 



{p n )i 



'So 



Si 



3V2. 



(18) 
(19) 



for the virial coefficients, since there is the dcutcron bound state in the 3 Si channel and the Levinson theorem requires 
<5c(0) = n c n with the number of bound states n c in channel c. Th e te rm — 3v2 from the lower boundary of the integral 
in b pn can be combined with the bound state contribution in ( 12 1 to give the modified single deuteron contribution 
(A d w ^A P An/2) 



-(d) 



3 / [J>p ~\~ Uf, 

exp 1 



A3 



T 



exp 



Bo 



T 



(20) 



Both expressions n[ d \ Eq. (14) and 
scattering states in the deuteron channel of 



(20) are related by a partial integration that changes the contribution of 
Eq. (15). Thus we conclude that the contribution of the scattering states 

Part of the continuum contribution can be 



can be reduced if the bound state contribution is redefined adequately 
transferred to the bound state contribution. Therefore the subdivision of the correlated part of the density into a 
bound state and a scattering part is ambiguous. This problem has been extensively discussed for Coulomb systems 
where for the bound state part the Brillouin-Planck-Larkin partition function has been introduced to avoid artificial 
singularities |19) . As a consequence it makes only sense to consider correlations in a particular channel defined by the 
corresponding quantum numbers without dividing into bound and scattering parts. Thus we can write 



with the contribution of the free nucleons 



(p) , (n) 

nfree = + "l 



and the contribution of the correlated nucleons 

„ _ n ( d ) 4. „„ _ n (PP) + „(««) + n {pn)T=0 , (pn)T=l 
"corr "1 ' 2 "corr ~ "corr ' "corr ' "corr 



(21) 
(22) 
(23) 



that can be split further into the four nucleon-nucleon channels considered above. This point of view corresponds to 
the physical picture where all many-body states are considered as correlations of nucleons. 
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C. Continuum correlations and quasiparticle shifts 



Coming back to the generalized Beth-Uhlenbeck formula where the single-particle contribution is given by the 
quasiparticlcs containing the self-energy shifts, the contribution of the continuum states ( 10 1 differs from the simple 
Beth-Uhlenbeck formula d5J|. The difference (expressed by the 2sin 2 ((5 c ) term) arises because part of the interaction 
is already accounted for in the quasiparticle shifts. As shown below in a model calculation, the Born approximation 
is fully accounted for by the Hartree-Fock (HF) self-energy shifts. Therefore, the introduction of the quasiparticle 
picture reduces the contribution of the continuum to the virial coefficient. On the other hand, the generalized Beth- 
Uhlenbeck formula describes the high density region as well where the bound state contributions are suppressed by 
Pauli blocking. 

The lowest order contribution (Born approximation) of the two-body scattering continuum to the EoS can be 
shifted to the contribution of free quasiparticles with Hartree-Fock mean-field energy shifts. In the physical picture 
with nucleons as the only degrees of freedom, we start from a non-relativistic Hamiltonian 



H = ^E{l)a\ ai + - V{l2,l'2')a{a\a 2 ,a v 



(24) 



where {1} denotes the momentum pi, spin a\ and isospin t\ characterizing the neutron or proton state. The kinetic 
energy is E(l) = p\/{2mx) and the potential energy contains the matrix element U(12, 1'2') of the nucleon-nucleon 
interaction. Sophisticated nucleon-nucleon potentials can be used that are fitted to available experimental data. For 
demonstration, we perform exploratory calculations based on the Yamaguchi separable interaction potential |20j . 



V(12,1'2') 



A 7 2 t 2 

fip 2 +7 2 {p') 2 +7 



,7V * 



T 2 ,T 2 / 



(25) 



with the relative momentum p = (p2 — Pi)/2 and the center-of-mass momentum P = pi + P2 (supposing m p — m n ). 
Values for the parameters A, 7 fitted to nuclear data are given in Ref. [2D] . 

According to the Beth-Uhlenbeck formula, the contribution of the nn scattering continuum to the density reads 
(c.f. Eqs. ([15} and dl6b) 



(nn) _ 



2 5/2 



exp 



2pn 

T 



dE 



exp 



(26) 



where E — p 2 /ra n is the energy of relative motion. For the Yamaguchi potential we have in the weak scattering limit 
(Born approximation) the phase shift 



T 



4. f 1 + ^ 2 



(27) 



With this expression we find 



3 r°° 



A3 "\TJ**TJ (1 + x 2 ) 



dx ivTn exp - 



(28) 



by introducing x — \J m n E/"/ 2 . 

On the other hand, we have from the free quasiparticle contribution in lowest order with respect to A (Hartree-Fock 
approximation) 



nT = 2 



d 3 pi 
(2ny 



^ ( p 2 /(2m») + A HF ( Pl )-^ ] . 



T 



with the self-energy shift 

A HF (pi)=^VW2,PlP 2 )| e 



exp 



( p\l(2m n ) + A HF (p 2 ) 



T 



(29) 



(30) 
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Without the exchange term, that can be neglected at low densities, and after expanding for small A HF , 
the nondegenerate case 



d 3 Pi 



exp 



p\/(2m n ) - u r 



T 



1 - 



T 



(n) 

— n\ ' 



(nn) 



we have in 



(31) 



with a term 



(nn) 



(nn) . 



in addition to the free neutron term 



» 



It is given by 



d 3 Pl 
(2tt)< 



exp 



pj/(2m n ) - u n \ f d 3 p 2 A 
T ) T J (2tt) 3 fl 



-i 2 



(P2-Pi) 2 /4 + 7 2 



exp 



p\l{2m n ) - fi n 
T 



Introducing relative and center-of-mass momenta we have with x = p/j 



(nn) 



2 exp 

2 3/2 



2^, 
T 



exp 



d 3 P 
(2tt) 3 

2m»A At 3 

T / 7T 2 T 



exp 



P 2 
4m„T 



d 3 p 
(2tt) 3 
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p2 _|_ ^,2 



exp 



TTln.T 



(32) 



(33) 



(1 + z 2 ) 2 



exp 



2 2 

7 z ar 



The expressions (28) and (33) coincide in the lowest order of the interaction strength A. The factor 2 [sin(5i5 )] 



occurring in the continuum contribution to the density in the generalized Beth-Uhlenbeck formula ( 10 1 produces a 
lowest order 0(\ 3 ) in the weak interaction limit. Thus, the lower-order contributions in the continuum part of the 
ordinary Beth-Uhlenbeck formula are transferred to the single quasiparticle term n^ u of the density. 

As an example, calculations have been performed for a Yamaguchi interaction fitted to nn scattering interaction 
(7 = 1.3943 fm -1 , A = 704.76 MeV fm 3 ) and for a weaker interaction with half the potential strength, i.e. A/2, see 
Fig. [T] with the total baryon number density as a function of the chemical potential /i = /!„ (neutron matter) in 
different approximations. The correlated part of the density, that is contained in the quasiparticle picture, almost 
completely reproduces the second virial coefficient in the case of lower interaction strength. For stronger interaction 
the quasiparticle shift accounts only for a part of the second virial coefficient. Note that the nn interaction is strong, 
and the di-neutron is almost bound. Therefore, a mean-field approach is not sufficient to account for the continuum 
contributions. 




FIG. 1: (Color online) Baryon number density n as a function of the chemical potential /i for T = 10 MeV in 
neutron matter. Model calculations for a Yamaguchi interaction with an effective range parameter 7 = 1.3943 fm -1 
are shown for two values of the coupling strength: A = Ao/2 (left panel) and A = Ao (right panel), where Ao = 704.76 
MeV. The ideal gas of nucleons (black dotted line) is compared with the correlated part obtained from the standard 

Beth-Uhlenbeck formula (red dashed line), the generalized Beth-Uhlenbeck formula (green dash-dotted line) 
containing the factor 2 [sin 6] , the difference in the density if quasiparticles are introduced (blue dash-dotted line) , 

and the total correction part of density (black solid line). 



Neutron matter contains no clusters because the interaction is not strong enough to form a bound state. Instead of 
the cluster virial expansion, we have the standard virial expansion where the second virial coefficient as a benchmark 
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can be directly related to the observed phase shifts. For details see [5][5T]. The introduction of quasiparticle states 
and the reduction of the scattering contributions according to the generalized Beth-Uhlenbeck approach [13] is an 
important issue to go to higher densities. 



IV. QUANTUM STATISTICAL APPROACH AND THE CLUSTER VIRIAL EXPANSION 



A systematic approach to the cluster expansion of thermodynamic properties is obtained from quantum statistics. 
The grand canonical thermodynamic potential 



J = -PQ = -TlnTr c -("-m")/t 



(34) 



where P is the pressure and f2 the volume, can be represented by diagrams within a perturbation expansion [55], see 
also [19] . We have 



/> = iTrM-Gf] ' fi ' ,X 



2Vt ./ A 



TrS A GA, 



(35) 



P = P 



1 f 1 dA 



20 Jo A 



or 




where A is a scaling factor substituting the interaction V by AV^. G\ is the free single-particle propagator that gives 



the ideal part of the pressure Pq. The full single-particle Green function G\ and the self-energy £>, are taken with 
the coupling constant A. Depending on the selected diagrams, different approximations can be found. In particular, 
the second virial coefficient for charged particle systems has been investigated, see Ref. [TJj]. [H] 

An alternative way to derive the equation of state is to start from the expression for the total nucleon density 



n Tl (T, n P ,iJ, n ) 



l 



duj 
2^ 



(36) 



where is the system volume, T\ = n, p, and summation over spin direction is collected in the factor 2. Both the 
Fermi distribution function and the spectral function depend on the temperature and the chemical potentials /^ p ,/x n 
not given explicitly. The spectral function S\(l,u>) of the single-particle Green function G\{l,iz v ) is related to the 
single-particle self-energy E(l,z) according to 



2ImEi(l,<j-zO) 



(lu - E(l) - ReSi(l, uj)) 2 + (Im£i(l, cu - iO)) 2 



(37) 



where the imaginary part has to be taken for a small negative imaginary part in the frequency. 

Both approaches are equivalent. As shown by Baym and Kadanoff [52], self-consistent approximations to the 
one-particle Green function can be given based on a functional $ so that 



£i(M') = 



6<f> 



(38) 



Different approximations for the generating functional $ are discussed in the App. [D] The self-consistent ^-derivable 
approximations not only lead to a fully-conserving transport theory. In the equilibrium case they also have the 
property that different methods to obtain the grand partition function such as integrating the expectation value of 
the potential energy with respect to the coupling constant A, or integrating the density n with respect to the chemical 
potential [A, lead to the same result. In particular, with 



also 



J = -Tr ln(-Gi) - TrE x Gi + $ 
ldJ 



(39) 



(40) 



holds in the considered approximation. 



10 



The latter approach using Eq. (36 1 has been extensively used in many-particle systems [13-15 , in particular in 
connection with the chemical picture. An analysis of the self-energy allows to work out a diagram technique that 
treats bound states on the same footing as "elementary" single particle described by the free propagator 

G[°\l,z) = 1— . (41) 

z - Eiipx) 



Considering the A-particle propagator ( A6 ) in the low-density limit where we can drop all medium effects, the solution 



of the Bethe-Salpeter equation in ladder approximation gives the propagator for the A-particle bound states 

G h A T d (l---A;l' ...A';z A ) = {l...A\^ AuP ) ^ ty AvP \V . . . A') (42) 

z A ~Ef v {P) 

where v indicates the internal quantum state of the A-particle bound state, after separation of the center-of-mass 
momentum P. As a new element, the bound state propagator is introduced as indicated in Fig. [2] This bound state 



propagator has the same analytical form like the single particle propagator (41 1, besides the appearence of the internal 
wave function that determines the vertex function. 



G A 



FIG. 2: Splitting of the A-particle cluster propagator into a bound and scattering contribution. Note that the 

internal quantum number has been dropped. 

As example, in App. [D] different approximations are obtained such as the nuclear statistical equilibrium (NSE) and 
the cluster mean-held (CMF) approximation using the chemical picture. These approximations are based on the bound 
state part of the A-particle propagator. They give leading contributions in the low-density, low temperature range 
where bound states dominate the composition of the many-particle system. From the point of view of the physical 
picture, these contributions arise in higher orders of the virial expansion of the equation of state. As example, the 
formation of the A-particle bound state is seen in the A-th virial coefficient, the mean-field shift due to a cluster B 
in the (A + B)-th virial coefficient. The chemical picture indicates which high-order virial coefficients of the virial 
expansion are essential, if the many-particle system is strongly correlated so that bound states are formed. 

As an example, Hydrogen molecules dominate the electron-proton system under normal conditions. The effective 
interaction can be calculated, and the corresponding virial coefficient determines the non-ideal part of the Hydrogen 
EoS. Local effective interactions are the Lennard- Jones or the Morse potential, and the Mayer cluster expansion can 
be used to evaluate the Hydrogen- Hydrogen virial coefficient. Within the physical picture based on electrons and 
protons, the eighth-order virial coefficient has to be analyzed to find the non-ideal features of the Hydrogen gas. 
Similar approaches can also be used for a matter in nuclear physics, introducing effective interactions such as the 
Ali-Bodmer potential. 

Using the quasiparticle approach, the EoS, Eq. ([8]), is obtained, see App. [Aj In contrast to the NSE, Eq. ([!]), 
medium dependent quasiparticle energies are used. In the cluster mean-field (CMF) approximation only the first- 
order terms of the cluster-cluster interaction are taken into account. The remaining part of the continuum correlations 
is neglected. 

To improve the approximation, the scattering part of the A-particle propagator has to be considered. It contributes 
also to the A-th virial coefficient. The scattering processes within the A-particle system can have different channels. 
As an example we discuss here binary elastic scattering processes between sub-clusters Ai and Ai of the system of 
A particles, A = Aj + A 2 . Binary phase shifts S AltA2 (E) are introduced that describe the corresponding scattering 
experiments. They can also be calculated within few-body theory. Besides the effective interaction between the sub- 
clusters that are depending on the internal wave function of the sub-clusters, also virtual transitions to excited states 
have to be taken into account. In general, the effective interaction is non-local in space and time, i.e. momentum and 
frequency dependent. 



The generalized cluster Beth-Uhlcnbcck formula ( 10 1 is obtained when in particle loops not the free propagator, 
but quasiparticle Green's functions are used. If the quasiparticle shift is calculated in Hartree-Fock approximation, 
the first order term of the interaction must be excluded from the ladder X , 2 addcr matrix to avoid double counting. The 
bound state part is not affected, it is determined b y an infinite number of diagrams. The scattering part is reduced 



subtracting the Born contribution as shown in Eq. (10) by the 2 [sin((S c )] 2 term; for the derivation see Ref. 



The continuum correlations that are not considered in the NSE give a contribution to the second virial coefficient 
in the chemical picture. We can extract from the continuum part two contributions: resonances that can be treated 
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like new particles in the law of mass action, and the quasiparticle shift of the different components contributing to 
the law of mass action. Both processes are expected to represent significant contributions of the continuum. After 
projecting out these effects, the residual contribution of the two-nucleon continuum is assumed to be reduced. One 
can try to parametrize the residual part, using the ambiguity in defining the bound state contribution, see Sec. |IIIB| 
Eventually the residual part of the continuum correlations can be neglected. 



V. CONCLUSIONS 



We have shown that it is possible to give an unified description for the nuclear matter equation of state that joins 
the region of saturation density, where quasiparticle approaches can be used, with the low-density region where the 
Nuclear Statistical Equilibrium or, more rigorously, the virial expansion can be applied. The chemical picture that 
allows for a systematic treatment of all bound states is used as guide line to formulate cluster expansions of different 
quantities such as self-energy, density and pressure. The inclusion of arbitrary nuclei as demanded by the chemical 
picture is indispensable to derive the nuclear matter EoS. 

As a main ingredient, the generalized Beth-Uhlenbeck formula [I2J E3 that treats the two-particle correlations 
already joins the low-density limit with the saturation density region. It is extended to clusters with arbitrary mass 
number |15j and has been investigated recently [4] to derive the thermodynamic properties of nuclear matter in the 
sub-saturation density region. 

Whereas the bound states are treated in a systematic way, the continuum of scattering states needs further in- 
vestigations, especially if many channels appear in the continuum. The approach presented is to extract important 
contributions from the continuum that are of relevance for the physical properties. We discussed three contributions: 

i) Resonances in the continuum can be considered similar to bound states. As example, we can consider A = 8 and 
the subdivision into two a particle clusters. The a — a elastic scattering phase shifts can be used to find a contribution 
to the corresponding cluster-virial coefficient. 8 Be as a sharp resonance can be treated as new particle in the chemical 
picture and should be projected out from the scattering contribution to the equation of state. 

ii) Mean-field effects can be extracted introducing quasiparticles not only for the single nucleon states, but also for 
the bound states. This reduces the interaction contribution to the continuum states. We have shown that in the 
weak scattering limit the continuum contributions are transferred to the quasiparticle energy shifts. In particular, the 
generalized Beth-Uhlenbeck formula shows this behavior. 

iii) There is no first principle distinction between the contribution of bound states and scattering states to the virial 
coefficient. For instance, partial integration and using the Levinson theorem gives different expressions for the contri- 
bution of bound states, see Sec. Ill B The contribution of the bound states to the virial coefficient is not unambiguously 
defined. This allows to find optimal approximations to take the contribution of the continuum into account. 

In conclusion, the cluster-virial expansion in combination with the consideration of excited states, resonances, and 
the introduction of the quasiparticle concept allows to extract relevant contributions from the continuum states. The 
remaining contributions of the continuum of the cluster-virial coefficients can be included into the sum over the internal 
quantum number v in the respective cluster-cluster channel. A systematic approach to the residual contribution 
that avoids double-counting is possible starting from the generating functional, but needs further study including 
investigations of topological aspects to characterize sets of relevant diagrams within the perturbation expansion. This 
is a complex problem and has to be worked out in future. Alternatively one can also use numerical methods to 
simulate the many-nucleon systems in the intermediate region. This way the interpolation between the saturation 
density and the low-density limit can be improved. 

Improved calculations of the nuclear matter EoS take the contribution of the continuum to the second virial 
coefficient into account. The ambiguity in the definition of the bound state contribution has been used in quantum 
statistical calculations, see the quantum statistical calculation in Ref. [3] and the recent work [5]. Formally, the 
summation over the internal quantum number v can be used to consider also the contribution of scattering states, 
as well-known from the Planck-Larkin partition function in plasma physics |19j . The effect of the correct treatment 
of continuum correlations for calculating the composition of nuclear matter is shown in Ref. [4]. The dissolution 
of clusters with increasing density was calculated within two different approaches, the quantum statistical approach 
where the second virial coefficient was taken into account, and the generalized RMF approach where continuum 
contributions arc neglected. The largest discrepancies are obtained for the deuteron fraction at high temperatures 
when the deuteron binding energy is small compared with the temperature (see also [5]). The influence of the correct 
treatment of continuum correlations in the deuteron channel on other thermodynamic quantities can also be seen when 
comparing both approaches. The systematic inclusion of further scattering phase shifts, e.g., in the a — a channel, 
within the cluster virial expansion would give an improved approach to all thermodynamic quantities. 
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Appendix A: Cluster expansion of the single-nucleon self-energy 

To include nuclei with arbitrary mass number A, a cluster decomposition of the single-nucleon self-energy can be 
performed [15 J. In particular, the inclusion of the light elements 3 H, 3 He, and 4 He has been discussed in [4]. We will 
not repeat the Green function approach here, see, e.g., Refs. [H[T5l[T6] but §^ ve om y some final results. The nucleon 
density is expressed in terms of the spectral function. The latter is related to the single-nucleon self-energy that is 
represented by Feynman diagrams. According to the chemical picture, the single nucleon propagators that occur in 
the self-energy have to replaced by propagators of arbitrary clusters. The cluster decomposition of the self-energy 
yields the law of mass action and the NSE. The mean-field approximation has to be replaced by a cluster mean-field 
approximation as given in App. [B] This way the quasiparticle concept for the single-nucleon state is extended to 
arbitrary clusters. 

The evaluation of the self-energy gives 

Z(l,z u ) = J2 Gj" ) _ 1) (2,..,i,« z ^^)T A (l...i I l'...A', Zyl ) (Al) 

A z A ,2...A 

with the free (A — 1) (quasi-) particle propagator 

G<°> u (2,...,A,z) = 1 fi,z 2 (2) f hZA (A) 

(A - 1)K ' ' ' J z-E 2 - ...-E A f A _ 1 . ZA _ 1 (E 2 + ... + E A ) 1 ; 

The quantity z v is the single-particle Matsubara frequency, and z A that of the A-particle system. The A-particle T 
matrix is obtained from the A-particle Green function by amputation. 

The A-particle Green function obeys in ladder approximation a Bethe-Salpeter equation (BSE) 

G A (1...A, I' ... A\ z A ) = (1...A, z A )5 lv . . . 8 AA , 

+ Y, Gf{l...A 1 z A )V A ' mi {l...Ay...A")G A {l"...A\l'...A' 1 z A ) (A3) 

l". ..A" 

where V AM (1...A, V ... A 1 ) = Y,i<j[Vij + AV^] is the interaction within the A -particle cluster, including mean-field 
contributions, see App. [5] The free A-quasiparticle Green function results as 



[1 - A ' ZA) ~ z A -ET(l) ET{A) ■ (A4) 



The approximation of an uncorrelated medium, see App. [Cj leads to the effective occupation numbers 

^"- oxp |£ Hi)/r->Vri + i "T(l?) 3,2 ^ r ' 1>/r < A5 > 

in the low-density, non-degenerate limit (jl T < 0), where jl T is determined by the normalization condition 2 ^2 p /i(p) = 
n r , where r denotes isospin (proton or neutron). 

The solution of the BSE is given by an expansion, bilinear in the wave functions 

G A (1...A, I'... A', z A ) = V i, AvP (l . . . A) 1 =^V^p(l' ■■■A'). (A6) 

^ z A E ■ - 



AuP 
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The summation over the internal quantum states v includes besides the bound states also the scattering states. The 
A-particle wave function and the corresponding eigenvalues follow from solving the in-medium Schrodinger equation 

[Er(l) + EZ(PMa„p(1 ■ ■ ■ , k, ■ ■ ■ A) 

+ E II S kk ^ AuP (l'...,k',...A')=0. (A7) 

V ...A' i<j k^i.j 

This equation contains the effects of the medium in the quasiparticle shift as well as in the Pauli blocking terms. 
Obviously the bound state wave functions and energy eigenvalues as well as the scattering phase shifts become 
dependent on the cm. momentum P, temperature T, and densities n p ,n„. 

Two effects have to be considered, the quasiparticle energy shift and the Pauli blocking. Detailed results have been 
obtained for the two-nucleon case. The shift of the binding energy and the medium modification of the scattering 



phase shifts are discussed extensively, see (T3l [40] . The solutions of the in- medium Schrodinger equation ( A7 ) for 
A = 2, 3, 4 was parametrized recently [IB] . 

Appendix B: The Cluster-mean field (CMF) approximation 

The chemical picture gives the motivation to extend the mean-field approximation for the case of cluster formation. 
Bound states are considered as new species, to be treated on the same level as free particles. A conserving mean-field 
approach can be formulated by specifying the Feynman diagrams that are taken into account when treating A-particle 
cluster propagation [18) . The corresponding A-particle cluster self-energy is treated to first order in the interaction 
with the single particles as well as with the f?-particle cluster states in the medium, but with full anti-symmetrization 
between both clusters A and B. We use the notation {A, v, P} for the particle number, internal quantum number 
(including proton number Z) and center of mass momentum for the cluster under consideration and {B, 9, P} for a 
cluster of the surrounding medium. 

For the A-particle problem, the effective wave equation reads 

[£7(1) + . . . E{A) - E AvP \ip AvP {l ...A) 

A 

V.. .A' i<j 

+ E y "m mf (1 ■■■A,l' ... A')^ AvP {l' ...A') = 0, (Bl) 
V.. .A' 

with Vtf (1 . . . A, V . . . A') = 7(12, 1'2% 3 / . . . 8 AA , . The effective potential 7^(1 ■ ■ ■ A, V ... A') describes the in- 
fluence of the nuclear medium on the cluster bound states and has the form 

K A m mf (1 .-.A, l'...A') = E A W^' • ■ ' 5aa ' + E'^ 1 ■■■A,l'...A'), (B2) 



with 



A(i) = E(^ 12 ' 12 )-/( 2 )-EEE E m e bvp)x 

2 B=2 v p 2...B V...B' 



x E • • • B, 1' . . . B')i> BP p(l . . . B)r BpP (l' ...&), 

i<3 



AV^(1 . . . A, 1' ... A') = - 1 + /(l'))V(12, 1'2') + (B3) 

oo B 

+ EEE E Ib(E BvP )Y j V 1 B j {12'...B',Y2...B)x 

B=2 vP 2...B2'...B' j 

x^ pP (22 . . . B)V Bp p(2'2' ...B>) U 33 , . . . S AA i , 
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/(l) = /i(l) + EEE fB{E B - uP )\^ B ,p{l ...B)\\ (B4) 

B=2 v p 2...B 

where (see Eq. ([2]); the charge number Z counts as internal quantum number) 

Ia{E) = e(E _ Zllp _ (A _ z)lln)/T _ ■ (B5) 

We note that within the mean-field approximation, the effective potential remains energy independent, i.e. 

instantaneous. The quantity /(l) describes the effective occupation of state 1 due to free and bound states, while 
exchange is included by the additional terms in AVj^ and A(l), thus accounting for antisymmetrization. 

Of course, the self-consistent solution of the cluster in a clustered medium is a rather involved problem which has 
not been solved until now. In particular, the composition of the medium has to be determined, with energy shifts of 
the different components (clusters of B nucleons) in the medium solving the effective wave equation for the i?-nucleon 
problem. 

Two effects have to be considered on the single nucleon level, the quasiparticle energy shift and the Pauli blocking. 
Phase space is also occupied by clusters as expressed by /(l). This effective occupation number is normalized to the 
total nucleon density. As approximation, a Fermi distribution function that is normalized correspondingly has been 
used in recent works [4]. Obviously the bound state wave functions and energy eigenvalues as well as the scattering 
phase shifts become dependent on temperature and density. 



Appendix C: Parametrization of single- nucleon quasiparticle shifts 

Different approaches to determine the single-nucleon quasiparticle self-energy shifts are known from the literature. 
We will not give an exhaustive review but mention only some general features, see also [2J. There are first principle 
approaches that start from realistic nucleon- nucleon interactions and solve the many-particle problem numerically. 
Quasiparticle self-energy shifts can also be obtained by different approaches such as the Dirac-Brueckner Hartree-Fock 
methods or phenomenological density functionals. The latter are based on effective density dependent interactions 
such as Gogny or Skyrme forces or use relativistic mean-field (RMF) concepts, see [U [25H3T] . Parameter are adjusted 
to reproduce nuclear bulk properties and properties of nuclei. 

There are several relativistic mean- field parametrizations actually used, like TW99 [3T], DD [35], DD2 pQ, TM1 
[7J TMA [5], FSUGold [33] and DDME<5 [35]. We focus on the density dependent DD2 parametrization g] 
adjusted to experimental properties of nuclei. In particular, this parametrization predicts reasonable values for the 
saturation density (n sat = 0.149 fm~ 3 ), binding energy (E/A = -16.02 MeV), compressibility (K = 242.7 MeV, what 
is not far from the experimentally determined value of 231 ± 5 MeV [35]), symmetry energy (J = 32.73 MeV, what 
fits the experimental value 31.3 MeV [37]), and the symmetry slope parameter (L = 57.94 MeV consistent with recent 
experimental constraints [38]). For neutron stars the model predict a maximum mass of 2.44 M Q , not in conflict with 
the largest known mass of M = 1.97 ± 0.04 M [59] . 

For direct use, a parametrization for the DD model [32] was presented in Rcf. We give here an improved 
parametrization of the DD2 model [1] in form of a Pade approximation. The variables are temperature T, baryon 
number density n = n n + n pi and the asymmetry parameter S = 1 — 2Y p with the total proton fraction Y p — n p /n. 
The intended relative accuracy in the parameter value range T < 20 MeV, n < 0.16 fm~ 3 is 0.001. 

The spectral function in the RMF approach gives the quasiparticle dispersion relation (i = p, n, no antiparticles) 
for the (non-relativistic) single quasiparticle energies 



e t (k) = /[m, - S(n, S, T)] 2 + fc 2 + V^n, 6, T) - m, . (CI) 
In the non-relativistic case we have 

e t (k) = -. If +Vi(n,5,T)-S(n,6,T) . (C2) 

2 [mj - S(n,5,T)\ 



With the quasiparticle energies, the chemical potentials fa follow from solving 

tt 2 J exp{[e 4 (fc) -^}/T} + l 



k 2 

dk — — — . (C3) 



The scalar self-energy (identical for neutrons and protons) is approximated as 

l + Sl {S,T)n + s 2 (6,T) n 2 



S(n, 8, T) = sq (5, T) n 



l+s 3 (S,T) n + s 4 (S 7 T) n 2 



with coefficients 



Si,o + Si.i T + s h2 S 2 + s ij3 S 4 
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(C4) 



(C5) 



densities n in fm 3 and temperatures T as well as the self-energies S, V in MeV. Parameter values are given in the 
Table D 







i = 


i = 1 


i = 2 


i = 3 


i = 4 


3 = 





4463.117 


20.56456 


15.98022 


24.27416 


114.5972 


3 = 


1 


-6.609841 


-0.040985 


0.866352 


-0.074176 


1.349746 


3 = 


2 


-0.170252 


-0.339370 


-2.020097 


-0.542662 


2.674353 


3 = 


3 


4.111559 


0.997156 


-3.018041 


1.196491 


0.726793 



TABLE I: Coefficients Sjj for the Pade approximation of the scalar self-energy S(n,6,T). 

The vector self-energy V p (n,S,T) = V n (n, —5,T) is approximated as 

V p {n, S, T) = v (S, T) n— — — — — - — 

1 + V3(d, Tj n + u 4 (d, T) n z 



with coefficients 



and parameter values given in Table [TTJ 



Vi = Vifi + ViA T + v it2 5 + v h3 S 2 



(C6) 



(C7) 





i = 


i = 1 


i = 2 


i = 3 


% = 4 


3=0 


3403.144 


0.662946 


10.77796 


3.432703 


23.01450 


3 = 1 


0.000052 


-0.006142 


0.004432 


0.000104 


-0.033018 


J — 2 


486.581687 


1.140795 


-0.802040 


1.548693 


5.922645 


i = 3 


-2.420361 


-0.717645 


0.457561 


-0.336038 


0.050892 



TABLE II: Coefficients Vij for the Pade approximation of the vector self-energy V p (n, S, T). 



Appendix D: Generating functional approach for the cluster virial expansion 



In this Appendix we present instructive examples of approximations to the $ functional on the level of Feynman 
diagrams to demonstrate the power of this approach. As shown by Baym and Kadanoff [21], any choice of a subset 
of diagrams for $ defines a selfconsistent approximation, where every Greens function is dressed by its appropriate 
self-energy obtained by the variation of $. The advantage of such so-called conserving or ^-derivable approximations 
is that they satisfy conservation laws and guarantee thermodynamic consistency. 

In the past the approach has been applied to a wide variety of many-body problems [23H28j but not yet to the 
problem of cluster formation in nuclear matter which we are going to discuss here. 

In the physical picture the generating functional is represented by the sum of all diagrams that consist of A = 2, 3, .... 
free fermion loops connected by arbitrary numbers of interaction lines. Exchange diagrams have to be added so that 
the correct symmetry is realized for identical particles. Disconnected diagrams have to be dropped. Only topologically 
different diagrams are allowed. With the self-energy that is obtained by opening one Green's function line, we can 
evaluate the density as done in this work. A cluster decomposition for $ is shown in Fig. [3] 

A general formalism in the chemical picture is not available at present. The full A-nucleon (cluster) propagator may 
be decomposed into a bound and scattering contribution ( 42 ) , given diagrammatically in Fig. [2] The concept is to 
consider the bound state part of the A- particle propagator on the same footing as the single-particle propagator. We 
illustrate this concept by considering special approximations that have been used in the present work. The relevant 
diagrams are selected by construction. 
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FIG. 3: Cluster decomposition of the generating functional $. Some contributions for A — 2,3 are shown explicitly, 
the ladder sum has to be completed considering an arbitrary number of interaction lines. The single particle loop 

has been added to represent the ideal part of the density. 



The ideal gas of nucleons follows from the first diagram of Fig. [3] The Hartree-Fock (HF) contribution to the 
second virial coefficient follows from the second diagram and the corresponding exchange term. The sum of all ladder 
diagrams with A — 2 yields a Bethe-Salpeter equation that contains a scattering part and eventually a bound state 
part. For A = 2, this binary approximation for <f> generates the standard Beth-Uhlenbeck formula. In the chemical 
picture, we obtain the nuclear statistical equilibrium (NSE) supplementing the single nucleon propagator (first diagram 
of Fig. [3]) by the propagator of bound states, see Fig. |lj 



n NSE : 

A 




FIG. 4: Generating functional for the nuclear statistical equilibrium (NSE). 

Within the single-nucleon approach, the corrections in lowest order of the interaction are given by the Hartree- 
Fock approximation. The corresponding Hartree-Fock approximation $ HF (as the simplest version of a quasiparticle 




FIG. 5: Hartree-Fock approximation for the generating functional $ qu . The dotted line is the interaction. The 
crosses lines denote antisymmetrization. The diamonds at the endpoints of Green functions denote amputation 

(multiplication with the inverse). 

approximation) is shown in Fig. [5] The dotted line represents the interaction. The crossed lines denote antisym- 
metrization. Note that bookkeeping has to be respected. Diagrams that are included in $ HF have to be subtracted 
from other groups. For instance, the first order diagrams of the ladder sum in Fig. [3] (the second and third diagram) 
have to be subtracted. 

Turning to larger clusters, we can also define a quasiparticle propagator according to Eq. ( |A6| . The presentation by 
diagrams using the A-cluster self-energy is given in Fig. [6j Then, $ HF is completed accordingto the chemical picture 
by the cluster mean-field approximation for the generating function $ CMF shown in Fig. km This approximation 
contains the cluster mean-field as described in App. [Bj 
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FIG. 6: Quasiparticle propagator for cluster with cluster mean-field self-energy. The crossed lines means the full 
antisymmetrization between the clusters A and B. The diamonds at the endpoints of Green functions denote 

amputation (multiplication with the inverse). 
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